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I. INTRODUCTION 


The vacuum instability in the supercritical Coulomb field is one of the important problem of quantum 
electrodynamics that is exhaustively studied in M- Main physical quantity related to this problem 
is the induced vacuum polarization 3, [Ml- New great interest to the vacuum instability in the su¬ 
percritical Coulomb field was revived in connection with the Coulomb impurity problem in graphene. 
The effective fine structure constant in graphene is large, which gives the new possibility to study the 
two-dimensional quantum electrodynamics in the strong-coupling (in fact, supercritical Coulomb) regime 
and the existence of charged Fermi quasiparticles in graphene makes experimentally feasible to observe 
the vacuum polarization in strong Coulomb field. 

From the physical point of view there are two (subcritical and supercritical) regimes, depending mainly 
on the magnitude of Coulomb field charge. Theoretically, charged impurity screening in graphene in terms 
of vacuum polarization were investigated in [isl - l^ and in comprehensive reviews . These studies 

have shown that the vacuum polarization charge density is localized at the potential center [UMlll in 
the subcritical Coulomb potential while the vacuum charge density induced by the supercritical Coulomb 
potential has the form c/r " [mil causing a modification of Coulomb law at large distances. This 
behavior could be expected on dimensional grounds: S(r) and c/r^ are the only dimensionally consistent 
possibilities due to the scaling invariance of the massless Dirac equation (in the absence of any intrinsic 
length scale) [11. It will be remembered that close to the so-called Dirac points, charged quasiparticle 
excitations in the potential of graphene lattice are massless Dirac-like fermions characterized by a linear 
dispersion relation psl - f^ and so a sing-le electron dynamics in graphene is described by a massless 
two-component Dirac equation [13, [H, llj“[lE [Hi HI] • 


New important results related to screening in graphene were obtained in |15ll20| with taking into account 
electron-electron interactions in a self-consistent renormalization group treatment. It turns out to be the 
system self-consistently rearranges itself so that electrons at large distances never feel a supercritical 
effective coupling and the subcritical (stable) situation is therefore protected. This conclusion agrees 
with expectations for the corresponding problem in the convenient quantum electrodynamics, where the 
vacuum polarization charge in super-heavy nuclei behaves in such a way as to reduce the supercritical 
charge of nucleus to the threshold value [lol (see also [3l| , where the problem was investigated for super¬ 
heavy nuclei in the presence of a superstrong constant uniform magnetic field). 

The wonderful quantum phenomenon was revealed in [H| : the induced current density in graphene 
in the field of a solenoid turns out to be a finite periodic function of the magnetic flux. The induced 
polarization current in the QED 2+1 with an Aharonov-Bohm potential for massive and massless charged 
fermions was studied in [H] . The induced electric current due to vacuum polarization in the AB potential 
was observed in [H | in “a quantum-tunneling system using two-dimensional ionic structures in a linear 
Paul trap”. 

The dynamics of charged fermions in the superposition of Coulomb and AB potentials is governed by 
a singular Dirac Hamiltonian that requires the supplementary definition in order for it to be treated as a 
self-adjoint quantum-mechanical operator. So we need to determine the self-adjoint Dirac Hamiltonians 
and then to construct the correct Green function of the Dirac equation in the superposition of Coulomb 
and AB potentials. In such a superposition the subcritical and supercritical regimes are determined by 
the magnitudes of parameters (as well as the relations between them) characterizing the Coulomb and 
AB potentials. 

A main feature of the supercritical and (at some magnitudes of parameters) subcritical regimes is a 
nonuniqueness of the self-adjoint Dirac Hamiltonian; there exists a one-parameter family of self-adjoint 
Dirac Hamiltonians specified by additional boundary conditions at the origin (35j | . This is a manifestation 
of a nontrivial physics inside the origin and an interpretation of self-adjoint extension parameters is a 
purely physical problem [H |. For example, in an AB field the magnetic flux within the interior of the 
vortex determines the effective Hamiltonian outside it; the extensions can be parameterized by nontrivial 
boundary conditions on the wave functions at the origin and different choices lead to inequivalent physical 
cases |36l |. We can determine the self-adjoint extension parameter in terms of the parameter R that is 
the finite radius of a real solenoid [H| . The self-adjoint extension method was used to determine bound 
states of massive fermions in the Aharonov-Bohm-like fields [H| and of a magnetic dipole moment in 
electric and magnetic fields generated by an infinitely long charged solenoid, carrying a magnetic field 


Here we study the vacuum polarization of charged massless fermions in the superposition of Coulomb 
and AB potentials in 2-1-1 dimensions. We calculate the induced charge density in the vacuum in the 
subcritical and supercritical regimes, for the first time, using the Green function of the two-dimensional 


3 


Dirac equation with Coulomb and AB potentials. The presence of AB potential in our model gives the 
possibility to estimate effects, which are due to the interaction of the electron spin magnetic moment and 
the Aharonov-Bohm magnetic field. Since the interaction potential is repulsive or attractive for different 
signs of spin projection this feature must be taken into account in the behavior of wave functions at the 
origin. 

We shall adopt the units where c = H = 1. 

II. INDUCED VACUUM CHARGE DENSITY 


A. Solutions and Green function of the Dirac Hamiltonian 


We remind that the Dirac y'^-matrix algebra is known to be represented in terms of the two-dimensional 
Pauli matrices 7 ° = < 73 , 7 ^ = isai, 7 ^ = ia 2 where the parameter s = ±1 can be introduced to label 

two types of fermions in accordance with the signature of the two-dimensional Dirac matrices [s^; for 
the case of massive fermions it can be applied to characterize two states of the fermion spin (spin ”up” 
and ’’down”) [i^. 

We also note that by Coulomb potential in 2-1-1 dimensions, we mean potential that decreases as 
l/r with the distance from the source, having in mind that in a physical situation (e.g., in graphene), 
although the electrons move in a plane, their Coulomb interaction with the external field of the pointlike 
charge of an impurity occurs in a physical (three-dimensional) space and the electric field strength of 
the impurity is a three-dimensional (not two-dimensional) vector. Therefore, the potential Ao(r) ~ l/r 
(and not Ao(r) ^ logr, as would be the case in 2 - 1-1 dimensions) does not satisfy the two-dimensional 
Poisson equation with a pointlike source at the origin. Besides, in real physical space, because of the 
existence of finite magnetic flux inside solenoid $ = 2ttB the singular term including the spin parameter 
appears in the form of an additional delta-function interaction of spin with magnetic field of solenoid 
H = (0, 0, il) = V X A = BTTS{r) in the Dirac equation squared. The additional potential —seBS{r)/r 
will be taken into account by boundary conditions. It will be noted that such kind of point interaction 
also appears in several Aharonov-Bohm-like problems [4l|, . 

The Dirac Hamiltonian for a fermion of the mass m and charge e = — eg < 0 in an Aharonov-Bohm 
Aq = 0, Ar = 0, Aip = B/r, r = _ arctan(y/a:) and Coulomb Ao{r) = a/enr, Ar = 0, 

= 0 , a > 0 potentials, is 

Hd = aiP2 - sa2Pi + 0-3171-eoAo{r), ( 1 ) 

where = —idfi — eA^^ is the generalized fermion momentum operator (a three-vector). The Hamiltonian 
dD) should be defined as a self-adjoint operator in the Hilbert space of square-integrable two-spinors 'l'(r). 
The total Dirac momentum operator J = —id/dLp + SO 3/2 commutes with Hd. Eigenfunctions of the 
Hamiltonian ([1]) are (see, 0 ) 

’ ( 2 ) 

where E is the fermion energy, I is the integer quantum number. The wave function T is an eigenfunction 
of the operator J with eigenvalue j = + s/2) in terms of the angular momentum I and 

hF{r) = EF{r), F{r) = ) , ( 3 ) 


where 


- . d I + + s /2 a 

n = ISO 2 ——h CTi- \- 03 m -, 

dr r r 


fj, = eoB. 


( 4 ) 


It will be noted that the massless fermions do not have spin degree of freedom in 2-I-I dimensions [dj, 
nevertheless, the Dirac Hamiltonian ([IJ keeps the introduced spin parameter. 

The induced current density due to to vacuum polarization is determined by the three-vector j)j(r), 
which is expressed in terms of the single-particle Green function of the Dirac equation as 


Jm(i’) 



( 5 ) 





4 


where C is the path in the complex plane of E enclosing all the singularities along the real axis E 
depending upon the choice of the Fermi level Ep. The Green function G can be expanded in eigenfunctions 
of the operator J. Since the induced charge density in the vacuum is divergent and thus needs the 
renormalization it is helpful first to consider the model with charged massive fermions. For such a model, 
the radial parts (the doublets) of above eigenvalues must satisfy the two-dimensional Dirac Muation ([3]). 
Then the radial partial Green’s function Gi{r,r'; E) is given by (just as in 3-1-1 dimensions [ 8 |) 

Gi{r,r'-E)-f° = :^^^[0(r'- r)C/fl(r)C/|(r')-h 0(r -/)t//(r)C/^(/)], ( 6 ) 

where W(E) is the (r-independent) Wronskian, defined by two doublets V and F as Wr(F, F) = Vi(J 2 F = 
(i'i /2 — fiV 2 ) and UR{r) and Ui{r) are the regular and irregular solutions of the radial Dirac equation 
{h — E)U{r) = 0; the regular (irregular) solutions are integrable at r —>• 0 (r —^ oo). We see that the 
problem is reduced to constructing the self-adjoint radial Hamiltonian h in the Hilbert space of doublets 
F{r) square-integrable on the half-line. 

Since the initial radial Dirac operator is not determined as an unique self-adjoint operator the additional 
specification of its domain, given with the real parameter ^ (the self-adjoint extension parameter) is 
required in terms of the self-adjoint boundary conditions. Any correct doublet F{r) of the Hilbert space 
must satisfy the self-adjoint boundary condition [H, [H, |4^ 

{F'‘{r)ia 2 F{r))\r=o = (/ 1/2 - / 2 /i)|r=o = 0. (7) 


Physically, the self-adjoint boundary conditions show that the probability current density is equal to zero 
at the origin. 

We shall apply as the solutions of the radial Dirac equation (|1]) the doublets found in 


Fr 


fR{r,j,E)\ 

gRlr,l,E)J^ ^ 


fi{r,l,E) \ 
gi{r,j,E) J ’ 


where 


fR{r,j,E) = 
9R{r,l,E) 


\/m + E 


X 

y/m — E 


{ARMaE/\+s/2,'fix) + CpMaE/X- 


{ARMaE/X+s/2,jix) — CpMaE/X 


-s/ 2 ,7(2;)) , 

-s/ 2 ,7(3^)) I 


Cr s^ — aE/X 
Ar V + ma/X ’ 


( 8 ) 


(9) 


//(r,7, A) = ^ {AlWaE/X+s/2nix) + ClWaE/X-s/2ni^)) ^ 

^ (Jj 

gi{r, 7 , E) = - {AiWaE/x+s/ 2 ,jix) - CiWaE/x-s/ 2 ,'Y{x)) , = {ma/X - svf . (10) 

Here 

x = 2Xr, X = x/m^-E^, 7 = - a?, v=\l + g + sl2\, (11) 

Ar, Ap Cr, Ci are numerical coefficients and the Whittaker functions Ma,b{x) and Wc,d{x) represent the 
regular and irregular solutions. 

For a? < "f is real, for a? > "f = iy/a^ — = ia is imaginary. The quantities q = x/v'^ — and 

Qc = R 7 = 0 are called the effective and critical charge, respectively; it is helpful also to determine 
Qu = - 1/4 7 = 1/2. 


B. Induced charge density in the subcritical range 

In the subcritical range for q < q^, 7 > 1/2, we can chose as the regular solutions only ones FR{r) 
vanishing at r = 0 ; for 0 < 7 < 1/2 (g„ < g < qc) the regular solutions UR{r) must satisfy the self-adjoint 
boundary condition 0 and should be chosen in the form of linear combination of the functions FR{r) 
and Fj{r) [H, 


C7fl(r) = FR{r) + ^Fi{r). 


( 12 ) 
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The (r-independent) Wronskian is easily calculated to be 


Wr(Ffi,F/) = W{E,j) = {gRfi - fngi) = -2 ArAi -yryTT , /x (13) 

1 (7 + 1/2 — s/2 — aEjX) u -\- raajX 

where r(z) is the Gamma function ( 4 ^ and, therefore, in the subcritical range the single-particle Green 
function is completely determined. One can show that the contribution into the renormalized induced 
charge density coming from range 0 < 7 < 1/2 is small for any therefore it is enough to consider the 
case ^ = 0 in the subcritical range. Thus, we should chose as the regular solutions the functions FR{r) 
for all 7 > 0 to obtain 




+1 +00 

E E 


S = —1 1 — — OC 


fifn + 9i9R 

27rsW(£;,7) ■ 


(14) 


After some calculations, we obtain 


trGj,(r, r'; E)j° — - E -i- G ^ ^ ^ -'^))^aE/\+i/2,-^{x)WaE/x+i/2,'i{x)+ 


1 — — C 


r (27 + i) 


-|-m^a/A(7 — aE/X)MaE/\-l/2,'yix)WaE/\-l/2,'fix) + Ex—{MaE/\+l/2,'yix)WaE/\+l/2,'fix)) 


(15) 


where now 7 = y/v'^ — a?, v = I -\- g -\- 1 / 2 . 

We note that the singularities of Gi,(r, r'; E) can be simple poles associated with the discrete spectrum 
(in the range —m < E < m), and two cuts (—00, —m] and [m, 00) associated with the continuum spectrum 
in the ranges \E\ > m [i^ . 

For the partial Green function in a Coulomb field in 3-1-1 dimensions, the path G may be deformed to 
run along the singularities on the real E axis as follows: G = G_ +Cp + G+, where G_ is the path along 
the negative real E axis (Rei? < 0) from —00 to 0 turning around at i? = 0 with positive orientation, 
Cp is a circle around the bound states singularities with —m < FI < 0 (if we chose Er = —m), and G+ 
is the path along the positive real E axis (Rei? > 0) from 00 to 0 but with negative orientation (i.e. 
clockwise path) turning around at FI = 0 Q. In the considered case in 2-1-1 dimensions the path G may 
be deformed in the similar way [s^ . 

One can show that the contour of integration G with respect to E can be deformed to coincide with 
the imaginary axis and we obtain: 


+ 00 

jo{r) = -e J ^trG^(r,r',iF;) 7 °. 


(16) 


Applying the following integral representation [47l| 

xT{2-f + l) 


MaE/\±l/2,'y{x)WaE/\±l/2,'yix) — 


r(l/2-h 7 - oF;/A T 1/2) 


00 

/■ 


—X cosh s 


[coth(s/ 2 )] 2 “^/^±iF 2 ^(x sinhs)ds, (17) 


it is convenient to represent the induced charge density in the form 


jo(r) =-^ ^ I dE I dte (2acos(2aFl/A) coth<F 2 -y( 2 Ar/sinht) — 

2Er 
sinht 


+ 00 


00 00 


’ 0 0 


suY{2aE/X)l! 2 ^{ 2 Xr/ smht) | , 


(18) 


where now A = y/'w? -|- E^, I^{z) is the modified Bessel function of the first kind and the prime (here and 
below) denotes the derivative of function with respect to argument. 

Let us write g, = [ij] + a = n-\- a, where [g] = n denotes the largest integer < g, and 1 > a > 0. Hence 
n = 0,1, 2,... for /r > 0 and n = —1, —2, —3,... for /r < 0. Since signs of e and B are fixed it is enough 

to consider the only case g > 0. Then, denoting = I + a + 1/2, i/_ = 1 — a -|- 1/2, 7 + = Jv\ — 
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7- = — a^, where here and in all formulas below I = I + n, we rewrite the induced charge density in 

the form 


JoW = I 


+ 00 


oo oo 


— 2 Ar coth t 


Z =0 


0 0 


(2acos(2aill/A) cotht(J 27 ^ (2Ar/ sinht) + l 2 -y_ (2Ar/ sinht)) — 

2Er 
sinht 


sm{2aE / X){l2^^{2Xr / sinht) + l2^_{2Xr/ sinht)) (19) 


emphasized the prime (here and below) denotes We note that jo is odd with respect to charge e. This 
expression is similar to the induced charge density in a pure Coulomb field obtained in [2^. It contains 
divergence and its renormalization should be carried out using the obvious physical requirement the total 
induced charge to vanish. Due to the nonzero mass m the renormalization can be performed by usual 
way in momentum space: 

^ OO OO OO 

jo(/3) = p(/3) = J drE^'^^joir) = ‘^'^ j dx j dt j dy^^^^e~^^°''^*-Joil3ysmht/2b)f{y,t), 

0 0 0 

f{y,t) = y sin(^t)(/2^^(y)+ /2^_(2/)) - 2acothtcos(Ait)(/2T,+ (y)+ / 27 _(y)). (20) 

Here /3 = \q\/m, x = E/m, b = vT+l?, y = 2bR/ sinht, R = mr, y = 2ax/b. 

Further calculations with this term for a < 1/2,a ~ 0 are similar to those described in detail in 
for the vacuum polarization in a pure Coulomb field in the subcritical range. We introduce 
the renormalized induced quantity in momentum representation as n{/3) = limA^oo[/o(/3) — Hm^a-s-o/o(/3)] 
with an ultraviolet cutoff jifl < A (see, [l^H^). Because the nonzero mass m is the only dimensionful 
parameter in the Green function the resulting dimensionless function n{f3) can depend only on the ratio 
/3 = q/m. Accordingly, it becomes just a constant in the massless limit m —>■ 0, which is denoted as 
Q = lim,ra^o ^(/3)- It is obvious that Q is the induced charge density localized in the point r = 0 in 
coordinate space. Therefore the induced charge density in coordinate space has the form Q = Q5{y). Let 
us calculate Q. 

As was shown in [d^ the induced charge density (as the series in terms of powers of a) for small a 
contains divergences only in the coefficients of the a and a/ terms. We give ni(/3) with the a term that 
reflects the linear one-loop polarization contribution: 


OO / OO OO 


^i(/3) = —X! [dti fdy f dx^^^^e ^™®‘'‘Jo(/3j/sinht/26)/i(y,t)-h acotht(e + e m) 

^ '=Oo \o 0 ' 

as well as the renormalized induced charge Qi in the first order of a: 

^ OO / OO 

/(it I / dysinhtln(l/ysinht)e-*'™"‘'‘x 

i=o { \{ 

X yt{dL+{y) +I 2 ^_{y)) - cotht(/ 2 i.+ (y) +l 2 ,^_{y)) + acotht{e~'^''+* + 

^ OO / oo 

= — ^ f dt i f dj/ln(l/y sinhi) 


1=0 


0 \o 


t sinh ^y (ye ^ *{l 2 ,^+{y) + 


--{t cosh te-y *)(/ 2 .+ (y) + l 2 u_ m 


acotht(e' 


-2yjf.t _j_ —2y—t' 


( 22 ) 


Integrating (l22)l . we obtain 

Qi = ^ ^(/ -b 1/2 -b 0 ) 1 /'(I + 1/2 + a) + {I + 1/2 - a)tp'{l + 1/2 - a) - 2 - ^?3) 

where ip{z) is the logarithmic derivative of Gamma function [d^. For a <C 1, we find 
Qi = eoTr/d -b ea7r(2 In2 -b 1 — 7r^/d)a^ r; ea7r(0.25 — O.Oda^). 


(24) 
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The first term in Eq. (IM)) coincides with result obtained in [H, [13, US] • We note that the contribution 
into Qi from AB potential arises in the presence of Coulomb field only, is small and has opposite sign 
compared with a pure Coulomb one. 

We have carried out long calculations and got the total exact induced charge in the subcritical range 
in the form 


Q — Qi + Qr, 


(25) 


where 


Qr = — Vim 

TT 


1^0 


ln(r( 7 + - za)r( 7 _ - ia)) + - ln(( 7 + - za)( 7 - - ia))- 


-(( 7 + - ia)^p{^+ - ia) + (j- - ia) 2 /’('y- - ia)) + ia- 


l + 1/2 


{I + 1 / 2)2 _ 

— ia{{l + 1/2 + oi)^!}'[I + 1/2 + o) + (1 + 1/2 — ct)ilj'{l + 1/2 — ct))] ■ 


(26) 


This expression at a = 0 is in agreement with result obtained in [ 2 ^; the coefficient of the term at 
a — 0 was also found in perturbation theory The induced charge Q determined by Eq. (ESI) is 

negative. 

It is worth to note that the vacuum charge density is induced by the homogeneous background magnetic 
field in the massive and massless QED 2 _|_i 


C. Induced charge density in the supercritical range 


In the supercritical range q > qd^ = ia) the stronger singularity of the Coulomb potential at the origin 
has to be regularized, therefore, we need to determine the self-adjoint Dirac Hamiltonians specified, for 
example, by self-adjoint boundary conditions ©• Then, we straightforward construct the Green function 
in the form ([B]) in which the regular solutions C//{(r), satisfying ([7]), have to be chosen in the form of linear 
combination of the functions Efj(r) and Fj{r). For this range (7 = ia), the above two solutions Efi(r) 
and Fi{r) become oscillatory with the imaginary exponent and it is convenient to use in this range the 
self-adjoint extension parameter 6 [H, , related to ^ by 

Ar _ 2 i 0 ( -b a{m -f E)/\ + isa _ r(2zcr) _ T{-2ia) _ 

^Aj \Fo J ly + a(m + F)/A — isa r(l/2 — s/2 — aF/A + ia) r(l/2 — s/2 — aE/X — zct) 

where tt > 6 * > 0 and a positive constant Eq gives an energy scale. 

The Green function has a discontinuity, which is solely associated with the appearance of its singularities 
situated on a second (unphysical) sheet ReE < 0,lm£' < 0 of the complex plane E at q > qc, these 
singularities are determined by complex roots of equation W{E, ia) = 0 and describe the infinite number 
of quasistationary (resonant) states with complex ’’energies” E = |E|e*'”. For massless fermions (m = 0) 
and O' <C 1 their energy spectrum was found in [d^ : 

Ek^g^s = ReE = Eq cos(r) exp(—fc/2cr + 9/a + 7 rcoth 7 ra/ 2 a), (28) 

where r Ri —(1 -|- s)/da -I- lmtjj{ia) + 7 r/ 2 . Eq. (I28|) contains an essential singularity. These quasi-localized 
resonances have negative energies, thus they are situated in the hole sector. For ct <C 1 the imaginary 
part ImE = tanrEk^g^s ReE is very small and, therefore, the resonances are practically stationary 
states [d^. For example, for a = 1/2, s = 1 r Ri (1 -I- 0 . 0 d) 7 r. 

Physically, the self-adjoint extension parameter can be interpreted in terms of the cutoff radius i? of a 
Goulomb potential. For this, for example, we can c omp are Eq. (1281) with the spectrum of supercritical 
resonances in the cutoff Coulomb potential [H, HO, [Hljand approximately derive 9 ^ cr[c(a) -f InFoi?], 
where c{a) does not depend on R. We note that the cutoff radius R rather relates to a renormalized 
critical coupling that is also characterized by a logarithmic singularity at mR <C 1 in massive case [HO, HO] 
The simplest way to include these resonances in the induced charge density is to carry out the integral 
in E from —00 to 0 along the path S taking into account the singularities on the second sheet. After some 
calculations, we represent the induced charge (electron) density (IS|) as the sum of contributions from the 






subcritical and supercritical ranges, which have to be treated separately 


dE 




e I dE ^ fi{r,'y,E)fu{r,'y,E) + gi{r,-f,E)guir,'y,E) 


2 Tri 


2 TTi 


C 


s= — l /= —c 


dE 

2 TTi 


^ E 

} ^ 




(.{fi{r,ia, E) + g‘j{r,ia, E)) 
sW{E, ia) 


sW{E,-f) 


= jsubir) +jsup{r). 


(29) 


For the supercritical range 7 = *ct, 0 > 0 > tt, the sum in second term jaup is taken over Z of > 
(Z + /i + s/2)^. Then the paths C, S can be deformed to coincide with the imaginary axis E. 

The hrst term in Eq. (1291) was calculated and explicitly represented in previous subsection. The second 
term is convergent and its contribution to the induced charge density can be directly evaluated at to = 0 . 
Having performed simple calculations we leads jaup to 


jsupir) — 


E 


sv 


s+1 


r dE 

37^2^2 ar(2ia)r(—2ia) J Euj{a) 


r(TO + (i-s)/2-TOF;/|F;|) x 


where 


xr( — 1(7 + (1 — s)/2 — iaE/\E\)WiaE/\E\+s/ 2 ,iai‘^\E\r)WiaE/\E\-s/ 2 ,iai‘^\E\r) 


2 ^e f‘2\E\Y^*‘' ’^ + iaE/\E\+isa T{2ia) T{-ia + {1 - s)/2 - iaE/\E\) 
\ / u iaE/\E\ — isaT{—2ia) T(ia{1 — s)/2 — iaE/\E\) 


(30) 


(31) 


Rewrite (2|if|/ifo)“^®‘^ as exp(— 2 i( 7 ln(|if|/ifo)). As far as the integrand (l30l) decreases exponentially 
at \E\ ^ 1/r and strongly oscillate at \E\ —>■ 0, the main contribution to the integral over E is given by 
the region |if| ~ 1/r. So in order to evaluate jaup we replace \E\ by 1/r in the log-periodic term of the 
integrand m and obtain 


Jsiip(^) — 


E 


si/®+^r(i(T -I- (1 — s)/2 -I- ia) 


87 r 2 p 2 ^ au!-(a)r(2ia)T(—2ia) 


T(—ia -I- (1 — s)/2 -I- ia) x 


tXJ 

/ dE 

— W_ia+a/2,^cr{‘^Er)W_,a-a/2,iai‘^^'^)^ 


(32) 


where 


CU- 


(cr) = 1 - c 2 ^e-^ 2 »<Tln(i;or) ^-»« + ts(T r(2TO) T{-ia + {1 - s)/2 + ia) 
V — ia — isa T{—2ia) r{ia -I- (1 — s )/2 -|- ia) 


(33) 


Because of the complex singularities on the unphysical sheet at g > (?c, the Green function and jaup(,r) 
are complex though for cr <C 1 their imaginary parts are small. In terms of the physics the complex Green 
function probably reflects the lack of stability of chosen (for constructing Green function) neutral vacuum 
for q > Qc (see, also Q). 

Now we can integrate in Eq. (15^ using formula [l^ 


OO 

J Ew_ia+a/2,iai‘^Er)W_ia-a/2,tai‘^Er) = 


s sin(27ri(T) 


1 


1 


r((l — s)/2 + ia + i(T)r((l -I- s)/2 -I- TO — ia) r((l — s)/2 + ia — TO)r((l -|- s)/2 -|- to -I- to)J^ ^ 
and after simple transformations we finally find the induced charge density in the supercritical range as 


J sup 


(r) = 


27r^r^ 


y Re- 
2 -^ / 


l,s:iy<a 


-(o')' 


(35) 


The main effect, arising at supercritical regime, is that the induced vacuum polarization for non¬ 
interacting massless fermions has a power law form (^ c/r^) whose coefficient is log-periodic func¬ 
tions with respect to the distance from the origin. In the subcritical regime the induced vacuum 
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charge is localized at origin and exhibits no long range tail. As an example, we consider Eq. (1551) 
for 1/2 — a < a < 3/2 + a, 1/2 a > 0, when just the lowest l + n,s channels are supercritical, and find 


fsupir) = 


27r^r 


crRe 


2 _ j^^g^2id+2itTln(Eor) 


(T — (T± 


1 - + A^[{a - tT)/(a + ’ 


(36) 


where 


^ r( 2 icr)r(—icr + ia) 

r(—2fCT)r(icr + ia) ’ 


Qi CT . “7 “7 

^ = 2-, h'± = 1/2 ± a^a± = a — 


(37) 


The induced charge density (I36|) resembles the local density of states which also exhibits resonances at 
the negative energies 

Using the known representation 


Argr(a: + iy) = y 


-c + E 



1 

y 



y 

cLL U L cLi. 1 . 

a; -I- n — 1 

)\ 


where C — 0.57721 is Euler’s constant, we finally obtain the induced charge density (1551) in the form 

■r ( . V R __ 

JsupV) 27r2r2 ^ ®1 - |A|ze2*®+2*‘’’M-Eoi-)+iV'+ |y4|2[(Q_ (jj/(a_l_crj]g4ie+4io-in(£or)+2iV'' ' 

Here 

, , . /^2cr „ 2cr 2ncr \ 

'll; = ArgA = —tt — 2Ccr + >-2 arctan-h arctan ^5 - 5 - . (39) 

^—' \ n n / 

n—1 ^ 

For small cr <C 1, Eq. (1551) takes the simplest form 

fsupir) = = v'a2-(l/2±a)A (40) 

Here “±” sign before a for hxed sign of a, in fact, corresponds to the spin projection sign. One sees that 
Jsup(^) i® odd with respect to the fermion charge e and even with respect to a. The contribution into 
the induced charge density due to the AB potential has opposite sign compared with a pure Coulomb 
one. It is of importance that the induced charge density BOl) at u ^ 1 does not contain at all the 

self-adjoint extension parameter 9. From the physical point of view, increased a near the point (7 = 0) 
the transition will occur from the subcritical range to the supercritical one, which can be symbolically 
characterized by 7 —^ cr and then a small change in a such that q > leads to a sudden change in the 
character of a physical phenomenon due to emerging of infinitely many resonances with negative energies. 
Put another way the character of a physical phenomenon itself must be due only to physical (but not 
mathematical) reasons. We also note that the expression jl^p{r) (at a = 0) is in agreement with results 
obtained in |15l| for the problem of vacuum polarization of supercritical impurities in graphene by means 
of scattering phase analysis. 


III. CONCLUDING REMARKS 

In this paper we have studied the vacuum polarization of charged massless fermions in Coulomb and 
AB potentials in 2-1-1 dimensions. In particular, we have calculated the induced charge density using the 
Green’s functions of the Dirac equation with the Coulomb and AB potentials. In subcritical regime the 
induced vacuum charge Q is localized at the origin and has a screening sign, leading to a decrease of the 
effective Coulomb charge; the contribution into Q due to the AB potential is small and has opposite sign 
compared to the Coulomb one. In the supercritical regime the induced vacuum charge, like the subcritical 
contributions, has a screening sign but it has a power law form, causing a modification of Coulombs law 
at large distances; the contribution into the induced vacuum charge due to the AB potential has opposite 
sign compared to the Coulomb one. 

Because a single electron dynamics in graphene is described by a massless two-component Dirac equa¬ 
tion we hope that our results can be applied in graphene with charged impurities. Furthermore, while 
the electron-electron interaction has been neglected, results of present paper may be useful to develop 
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further insight into the screening of the Coulomb impurity with taking into consideration of the electron 
spin and electron-electron interaction. To approach this problem one can write the self-consistent renor¬ 
malization group equations in the Hartree approximation in the subcritical range in the same spirit as in 
and within the Thomas-Fermi method for the supercritical range. We shall defer the self-consistent 
renormalization group analysis to a future work. 
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